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1. Introduction 

After many years of measurements, begun by H. Cavendish, the Newtonian gravi- 
tational constant value, said G , is still affected by a large error Gl. The CODATA 
recommends G = (6.6742 ± 0.001) x 10~^^^jt, meaning a relative standard uncer- 

tainty of 150 parts per million (in the following ppm) . This value has had recent 
confirmation by means, on one hand, of a super-conducting gravimeter ^ and, on 
the other hand, of a careful analysis of the possible beam balance nonlinearity SI. 
The original Cavendish method of measure, employing the torsion balance, still re- 
veals a large discrepancy from the recommended value: about 500 ppm IS! . This is 
probably due to imperfections of the crystalline structure of the torsion fibre l ^ l ^ l ^ i 
Moreover recent studies point out that the period of a torsion pendulum might vary 
under disturbances of environmental noise factors, see see Other authors suggest 
a possible deviation from Newton's law specified as an additional contribution of 
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Torsion Balance 
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Fig. 1. Schematic view of the torsional balance 



Yukawa potential type This paper first analyses the linear and non linear equa- 
tions governing the torsional balance rotational dynamics (Section [21). By means 
of a gravitational torsion balance, same values of G are obtained and summarised 
(Section [3]). The oscillation period's variation, due to a non linear effect, is then 
discussed (Section |4]). 



2. The basic equations 

The form of Newton's law of gravitation is 

F.alf. . (1) 

where G is the gravitational constant , M the great mass , m the small mass and 
r their relative distance. The Leybold balance represents a widespread instrument 
to determine the constant G, see Figure [TJ and is constructed with the following 
components. 

(1) A freely oscillating horizontal bar , of length 2d, holding two small lead balls of 
mass m as in Figure [1] supported by a torsion fibre that has a torsional constant 
r . 

(2) Two larger balls of mass M that can be positioned next to the small balls as in 
Figure [TJ The center of mass of the two ni and M are supposed to be all on a 
plane perpendicular to the fibre. 
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Fig. 2. Top view of the Cavendish balance 



(3) A luminous source directed toward the center of mass of the bar where is re- 
flected by a mirror. 

(4) A scale at distance / where the reflected light beam is measured. 

Thus the equilibrium position about which the pendulum oscillates is different for 
the two positions and it is this difference which we use to determine G. The Figure [2] 
reports a plot of the motion. The moment of inertia of the bar,/ is 

/ ~ 2m(f . (2) 

The fundamental equation of rotational dynamics is 

I^=Mg + M, + Mt , (3) 

where 

M, = -P^ , (4) 

Mt = -T^ , (5) 

here (3 is the coefficient of viscosity of air , Lp the angle between bar and bar itself 
when the torque is zero. This angle is measured in the anti clockwise direction. The 
term Mg represents the torque of the gravitational forces. From equation ([1]) , we 
obtain 



Mg = 2dGF{ip) 



(6) 
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where is a function of the angle (p. The law of dependence of F with ip is complex 
and will here be analysed. When the motion starts the resulting force is 

F(^) = FiM-i^2M , (7) 

where 

, , cos 03 , , 

Fi((^)=mM^ , (8) 

^/ N , rCOsff - ((y5 + a)) , sin((z) + a) 

F2 {ip) = mM )^ = mM — , (9) 

now a :— arcsin(^^-^) (where res are defined as in Figure[T]). With our data , see 
Section [3] the maximum angular excursion of the angle p is 

= arcsm( — ) . (10) 

The angle (p has a low values, see Table [1] and a Taylor series expansion that keep 
terms to order ip'^ will be adopted. This means to forget quantities less than 3 10^^. 
The series representation gives 

Ml = ^/{r + s)'^ + {d- dcosp)'^ ''(1 + '^(y + + ^'^^))^ ■ (H) 
Developing the last term with a Maclaurin series we obtain 

m^ril-^i^ + ^p)-f^^'+Oiip'))^r + d^ . (12) 
Z r r Ir 

As a consequence 



(r + dp)'^ 

and expanding the denominator we obtain 



1 _ 

Fi - mM- , (13) 



{r + dp)-^ = {r-^-2-p + i-p' + 0{p^)) , (14) 



-2 _ /„-2 o " ,„ I Q _ ,„2 

that means 



mM, 2d ,3^2 1 

An expression for can be obtained from equation ([9]) 

sin cos a + cos sin a r + d.^-|(p2 
1^2 ^mM ^ ^ mM ^ 3- . (16) 

u {d^ + r"^ ~2rdp - d^p"^)^ 

On Taylor expanding the denominator 

iAd-+r^-2rdp-d\r- = i^d'+^'rHl+J^^+'^^f^^'+Oip-)) 

(17) 
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and therefore 

(4^2 + r2)2 

Now _F from equation ([7]) can be expressed hke a Taylor expansion truncated at 



- + + ^2^^ , (19) 



where 



^2 = 3 1/2^ + 1/2 ^ — >- . 22 

' r2 ' (4^2 4_ ^2)7/2 

Three methods that allow to obtain an expression for G in terms of measurable 
quantities are now introduced. Further on the well known formula for G extracted 
from the Leybold manual is reviewed. 



2.1. Averaged G 

Let F{lp) — F for all the experience ; in first approximation we may assume that F 
is given by equation (|19p to first order 

F = A„ + A{^ , (23) 

where ^ is the average of the values that assumes between the first position , Pi , 
and the last position ,Poo, of the balance ; Ao,Ai are given by equations ( [20)1 and 
(|2T|) . The differential equation that describes the motion is 



Iip+f]f +Tip = 2dGF , (24) 

and it's solution is 

ipit) = ce-^* cos(wt + (/)) H , (25) 



where c represents the amplitude and 



*:=4 , (20) 



-^^^ . (27) 
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The angle (pac, that represents the bar position at the end of the phenomena Poo 
can be determined as follows 

2^ = ^^ ' 

and should be the same as predicted by the theory 

hm ip{t)^ ; (29) 

t—>-\-oo T 

therefore 

. (30) 
2dF ^ ' 

In order to continue a value for r should be derived. This can be obtained from the 

period of oscillation of the bar 

^ _ 27r _ 47r/ 

" w ^ - /?2 

We continue by identifying T with the empirical value T. We continue on assuming 
that ^ is small ; therefore from equations ([2]), ([50)1 and the following is 

obtained 

G = . (32) 

d(Ao + Ai^)r 



2.2. G wifft. air viscosity 

From formula (|4]) is possible to deduce the viscosity of the air once the coefficient 
of damping d is known, see Section [3] on data analysis. From formula ([30)1 and ((26|) 
we should add to the value of G reported in equation ((32| 



obtaining 



8IdF ^ ' 



d{Ao + Aiip)T^ 8Id{Ao + Aiif) 



^" + . (34) 



2.3. G to f/ie yirsf order 

Let assume that -?"((/?) is not constant, we can assume at the order 0{ip'^) with the 
aid of formula 

F^A^ + A^ip , (35) 

where Aq e Ai are defined in equations (PO)) and (PT|) respectively. In this case the 
law of motion is still equation ([3|) 

/ (/3 +(t - 2dGAi)v? = 2dGAo , (36) 
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and the solution is 



ip{t) = ce~** cos{uj't + </>) + 



2dGAo 



T - 2dGAi ' 

where the angular velocity uj' has now the following expression 



As a consequence 



' - 21 



T = ^ + 2dGAi + , 
4/ 



(37) 



(38) 



(39) 



and therefore 



Once equation 



2dGAo 



G = 



2d{Ao + Aiip^) 
is substituted in this relationship we obtain 



G = 



27r^/y?oc 
dAof^ 



SdlAo 



(40) 



(41) 



(42) 



2.4. G from Leybold manual 

The deduction of G through the Leybold torsional balance is widely known , see 
We simply report the final expression 



G=^g^x(l+/3) , 



where 



(43) 



(6^ + 4rf2)^62 ^4d2 



The meaning of the symbols is 

• 6:Distance between center of the great mass and small mass 

• A S:Total deflection of the light spot 

• d:The length of the lever arm 

• ZiDistance between mirror and screen 

• AI : Great mass 

• T Period of the oscillations 



(44) 
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parameter 


value 


unit 


M 


(1.5 ± 10-3) 


Kg 


m 


(1.5-10-2 ± 10-3) 


Kg 


r 


(4.65-10-2 ± 10-3) 


m 


d 


(5.0-10-2 ± 10-3) 


m 


I 


(5.475 ± 10-3) 


m 


V> 


(-6.715-10-3 ± 1.3-10-5) 


rad 


P 


( 1.432-10-'^ ± 1.11-10-*) 


s 



parameter 


value 


unit 


Ao 


(8.188-10-2 ±1.11-10-"*) 


m 


Ai 


( 0.1470 ± 2.85-10-4) 


m 


T 


(552.98 ± 0.16) 


s 


T 


(1047.0 ± 4.8) 


s 



method equation value uncertainty [ppm] accuracy %] 

G averaged §21 (6.67 ± 0.34)-10-"^ 52433 0.0161 

G with air viscosity jMl) (6.72 ± 0.35)-10-"^ 52433 0.69 

G to the first order gj) (6.80 ± 0.34)-10-"^ 49989 1.92 

G from Leybold manual (gS]) (6.71 ± 0.33) -10-"t^ 49600 0.64 



3. Analysis of the data 

The physical parameters as well their uncertainties are reported in Table [H The 
data were analysed through the following fitting function 

y{t) = Ao + Aicos{^)exp{--) . (45) 

The data has been processed through the Levenberg-Marquardt method ( subrou- 
tine MRQMIN in^l^ in order to find the parameters Ao,Ai ,T and r. The results 
are reported in Table [2] together with the derived quantities. 

The value of G can be derived coupling the basic parameters of the torsion 
balance , see Table [H and the measured parameters of the damped oscillations , see 
Table [21 Table [3] reports the four values of G here considered with the uncertainties 
expressed in absolute value and in ppm ; the precision of the measure in respect 
of the so called "true" value is also reported. A considerable source of error is the 
uncertainty in the determination of the span between t he t wo spheres that in our 
case is « 10-3m . Adopting a rotating gauge method the uncertainty in the 
determination of the span between the two spheres is w 0.5 x IQ-^m; this is the 
way to lower the uncertainty in Table [3l 
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4. Non linear effects in the vacuum 

By starting from the equation of rotational dynamics up to the second order 

lip+p(p+{T-2dGAi)ip-2dGAo = 2dGA2ip^ , (46) 

the case of /3=0 is analysed , 

I^+{T-2dGAi)ip-2dGAo^2dGA2ip^ , (47) 

that corresponds to perform the experiment in the vacuum. On dropping the con- 
stant term and dividing by / we obtain 

.. (T-2dGAi) 2dGA2 2 

J i^^=— _ (48) 

On imposing 

2 {T-2dGAi) 



(49) 

the nonlinear ordinary differential equation , in the following ODE , has the form 

, 2dGA2 . , , 

On adopting the transformation T ^ t ujq the nonlinear ODE is 

V +(p - eifi'^ = , (51) 

where 

_ 2dGA2 

' " (r - 2dGA^) ■ ^^^> 
The solution of equation ()51|) is reported in the Appendix [Appendix A| and in our 
case e=0.0187. We now have a measured period ,Tms, that is equahsed to the non 
linear value , Tjvl . The period of the linear case , , can be written as 

Tatl = Tm5 = l-OOOMTi , (53) 

and therefore 

1.00014 ^ ' 

In the various formulae of G without damping , for example equations p2p and (I43p , 
the periods Tj, and Tms are raised to the square 

" 1.00029 ' ^^^> 
and in the denominator, making the non linear Gnl greater than the linear Gl 

GjvL = 1.00029 Gl . (56) 

The value of this correction , SG, can be evaluated as a difference between 1 and 
the multiplicative factor of Gl 
6G= (1.00029- l)x 6.6742 Xl0-"m3fc5~\s"2 = 0.19 x 10"^^ m^fcg-^s-^ . (57) 

The official error on G is 0.1 x 10^^'^ m?kg^^s^^ and therefore the nonlinear cor- 
rection can be expressed as the double of the official error on G. 



2, 2008 16:24 WSPC/INSTRUCTION FILE cavendish 



10 Rossi Zaninetti 

Appendix A. The eardrum equation 

The equation 

x + x + , (A.l) 

is weh known under the name "eardrum equation". It can be solved , see^^, trans- 
forming it in 

n'-^XiT) + X{T) + e{XiT)f^O , (A.2) 
and adopting the method of Poisson that imposes the foUowing sohition to X 

X (T) = xo (T) + XI (T) e + X2 (T) , (A.3) 

and to r2 

VL = l+ uJie + oj2f? ■ (A. 4) 

The computer algebra system (CAS) gives 

wi = L02 := -5/12 . (A.5) 
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